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Abstract. We examine the generation of gravity waves by the surface convection zone of low-mass main sequence stars with 
solar metallicity. It is found that the total momentum luminosity in waves rises with stellar mass, up to the quasi-disappearance 
of the convection zone around 6500 K (corresponding to a mass of ~1.4 M Q for solar metallicity) where the luminosity dras- 
tically drops. We calculate the net momentum extraction associated with these waves and explain how the calculated mass 
dependence helps resolve the enigma of the Li dip in terms of rotational mixing, forming a coherent picture of mixing in all 
main sequence stars. 
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1. Clues to angular momentum transport in 
low-mass stars 

In many locations of the Hertzsprung-Russell diagram, stars 
exhibit signatures of processes that require challenging model- 
ing beyond the standard stellar theory 1 . In this context, rotation 
has become a major ingredient of modern models, especially 
when abundance anomalies have to be accounted for. In or- 
der to correctly describe the effects of rotation on stellar struc- 
ture and evolution as well as on their "byproducts" like surface 
abundances or chemical yields, special emphasis has to be put 
on the evolution of the angular momentum distribution; this in- 
deed is the feature that determines the extent and magnitude of 
rotation-induced mixing in stars. 

In recent theoretical developments, the internal rotation law 
evolves as a result of contraction, expansion, meridional cir- 
culation, shear turbulence and mass loss; mixing of chemicals 
is directly linked to the rotation profile (see e.g Talon 2003). 
The most sophisticated treatment of these hydrodynamical pro- 
cesses is based mainly on the work by Zahn (1992), Maeder 
(1995), Talon & Zahn (1997) and Maeder & Zahn (1998). It 
rests on only one assumption which is that the turbulence sus- 
tained by the shear is highly anisotropic. Two free parameters 
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1 By standard we refer to the modeling of non-rotating, non- 
magnetic stars, in which convection is the only large-scale mixing 
considered. 



describe the magnitude of the horizontal shear and the erosion 
of the restoring force due to both the thermal and the mean 
molecular weight stratifications. 

Such a self-consistent treatment has been successfully ap- 
plied in various parts of the HR diagram. The strength of the 
theory is that the use of the same free parameters give very 
satisfactory results for stars over a large range of masses and 
evolutionary phases. For example, it explains the observed He 
and N enrichment in main sequence O-type and early B-type 
stars, in OB supergiants (e.g., Gies & Lambert 1992, Lennon 
et al. 1991, Herrero et al. 1999) and in A-type supergiants in 
the SMC (Venn 1999), as well as the B depletion in main se- 
quence B-type stars (e.g., Venn et al. 1996); it also helps repro- 
ducing the number ratio of blue to red supergiants in the SMC 
(e.g., Meylan & Maeder 1982) and the observed WR/O ratios 
at solar metallicity. We refer to Maeder & Meynet (2000) for a 
more detailed description of the effects of rotation in massive 
stars and for relevant references. 

In low-mass stars, the strongest observational constraints 
come from the so-called Li dip and from helioseismology. The 
Li dip is a characteristic feature of lithium abundances which is 
seen both in the field and in open clusters (e.g., Wallerstein et 
al. 1965, Boesgaard & Tripicco 1986, Balachandran 1995). It 
corresponds to a narrow region in effective temperature (around 
T e ff ~ 6650 K) where the surface Li abundances are reduced by 
up to ~2.5 dex. In the same region, low-mass stars are begin- 
ning to be spun down efficiently early on the main sequence 
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(Boesgaard 1987), most probably via magnetic torquing. In 
modern rotating stellar models in which the transport of an- 
gular momentum and of chemicals by meridional circulation 
and shear turbulence is self-consistently taken into account, 
this leads to important lithium (and beryllium) destruction at 
the correct r eff (Talon & Charbonnel 1998, hereafter TC98; 
Palacios et al. 2003, hereafter PTCF03). In addition rotation 
induced mixing inhibits the atomic diffusion in a way that ex- 
plains the constancy of the CNO abundances within the Li dip 
(e.g., Varenne & Monier 1998, Takeda et al. 1998). Last but not 
least, these models do explain the evolution of Li surface abun- 
dances in evolved stars that originate from the hot side of the 
dip. It is worth recalling that all these results for low-mass stars 
on the blue side of the dip are obtained with the same values for 
the two above-mentioned parameters as those used in massive 
stellar models. 

For still cooler main sequence stars however, the magnetic 
torque strengthens as the stellar convective envelope grows. If 
we assume here that all the momentum transport is assured by 
the wind-driven meridional circulation as it is on the blue side 
of the Li-dip, then too much lithium burning is obtained be- 
low r e ff ~ 6550 K. Thus the rise of the Li abundance on the 
red side of the dip can be interpreted as the signature of an- 
other mechanism that efficiently transports angular momentum 
(but not chemicals) in lower mass stars (TC98). As a result, in 
these objects the magnitude of both meridional circulation and 
shear turbulence is reduced, as well as the Li depletion due to 
rotational mixing. 

Such a mechanism is also required to shape the Sun's flat 
rotation profile (Chaboyer et al. 1995, Matias & Zahn 1998). At 
the solar age, models relying solely on turbulence and merid- 
ional circulation for momentum transport still predict large an- 
gular velocity gradients which are not present in the Sun. 

Only two mechanisms have been proposed so far to trans- 
port angular momentum in addition to the classical hydrody- 
namical processes in order to enforce the flat rotation profile 
measured by helioseismology : magnetic fields (Charbonneau 
& MacGregor 1993, Barnes et al. 1999) and gravity waves 
(Schatzman 1993, Kumar & Quataert 1997, Zahn, Talon & 
Matias 1997). As explained in TC98, in order to obtain a con- 
sistent picture of rotational mixing in all types of stars, the cor- 
rect mechanism must become efficient in the center of the Li 
dip, namely at an effective temperature of ~ 6650 K. 

While the original propositions for momentum extraction 
by gravity waves have received much criticism (Gough & 
Mclntyre 1998, Ringot 1998), recent calculations performed 
by Talon, Kumar & Zahn (2002, hereafter TKZ02) show that, 
through differential filtering, gravity waves are indeed able to 
extract momentum from the solar interior. Here we explore the 
same physics in the stellar mass range around the Li dip in 
order to verify whether it has the proper mass dependence as 
described by TC98 (see Fig. 1). Our goal is to present a coher- 
ent understanding of rotational mixing in all stars of the main 
sequence. We will firstly recall the main descriptions that have 
been made for wave excitation (§ 2) and explain how thermal 
dissipation of these waves leads to the formation of a filtering 
shear layer (§3). Results for momentum extraction by waves in 
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Fig. 1. Measured Li abundances in the Hyades. Superposed is 
the pattern of efficiency of momentum transport required to 
produce the Li gap via rotational mixing 

stars of various masses are then presented (§ 4) and discussed 
(§ 5). 

2. Turbulent wave excitation 

2.1. Two excitation mechanisms 

The generation of gravity waves depends on the structure of 
the stellar convective envelope. On the main sequence, this re- 
gion is rather large in solar mass stars (~30% in radius), while 
it shrinks as the stellar mass rises, representing only ~1% in a 
1.45 M Q star lying inside the Li dip. The characteristic convec- 
tive frequency varies from v c ~0.75 /iHz to v c ~600 //Hz while 
the spherical harmonic number I corresponding to the largest 
convective scale varies from ( c ~60 to ( c ~1000 in the same 
mass range. These variations are not smooth, rapid changes of 
convective properties occurring when M, >1.3 M . We wish to 
investigate how the wave spectrum produced by these convec- 
tion zones varies as a function of stellar mass. 

The exact properties of wave spectra remain somewhat un- 
certain. Excitation can be related to internal stresses that corre- 
late with the mode's eigenfunction (Goldreich & Kumar 1990, 
Balmforth 1992, Goldreich et al. 1994). Gravity waves may 
also be excited by penetration below a convection zone, as ob- 
served in laboratory experiments (Townsend 1958) and in nu- 
merical simulations (Hurlburt et al. 1986, 1994, Nordlund et 
al. 1996, Kiraga et al. 2000). In a real star, both sources would 
contribute to wave generation and they are thus additive. 

In this paper, we examine both sources, with emphasis not 
on their absolute values, but rather on their variation with stel- 
lar mass. Our computations are performed for stars between 
1.1 and 1.45 M Q with Z = 0.02 (see Table 1). We use the same 
code and input physics (eos, opacities, nuclear reactions) as in 
PTCF03. Wave characteristics are computed on the basis of 
ZAMS stellar structures. Note that in these objects the depth 
and structure of the convective envelope does not vary signifi- 
cantly over a main sequence lifetime (at least up to the age of 
the Hyades). 
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Table 1. Characteristics of the stellar models (with Z = 0.02) on the ZAMS and maximum frequency v max used to capture the 
wave - mean flow interaction. l c is the largest convective scale, v c , the characteristic convective frequency and Kj, the thermal 
diffusivity. 
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2.2. Excitation by the Reynolds stresses 

First, we consider excitation by the Reynolds stresses as de- 
scribed by Goldreich et al. (1994). This description, which 
has also been used by Kumar, Talon & Zahn (1999, hereafter 
KTZ99), uses a free parameter that has been calibrated on the 
solar p-mode spectrum. 

The energy flux per unit frequency due to Reynolds stresses 
Te R is then 

mpl-^IM^ ^^ , (1) 

where f r and [((( + l)] 1 '^ are the radial and horizontal dis- 
placement wave-functions which are normalized to unit energy 
flux just below the convection zone, v is the mixing length con- 
vective velocity, L is the radial size of an energy bearing turbu- 
lent eddy, ti x L/v is the characteristic convective time, and 
h w is the radial size of the largest eddy at r with characteris- 
tic frequency of to or greater (h^ = Lmin{l, (2wtz,)~ 3 ' 2 }). The 
gravity waves are evanescent in the convection zone, the re- 
gion where they are excited. Their wave-functions <f,. and £/, are 
thus proportional to k~^ 2 exp(i J drk^. The above equation 
was derived under the assumption that the turbulence spectrum 
is Kolmogorov. The momentum flux per unit frequency Tj R is 
then 

f/(m,Aw)=-f/Kw). (2) 
to 

The momentum spectra corresponding to the various stellar 
masses are shown in Fig. 2 where Xy = Anr^^Fj is the momen- 
tum luminosity; note that some of the lower frequency waves 
have zero magnitude. This is due to the fact that the correspond- 
ing damping is too large to permit them to be waves. 

Spectrum characteristics evolve together with the structure 
of the convection zone which depends on the stellar mass. Low 
frequency waves disappear when the stellar mass increases. 
This is related to the fact that the surface convection zone then 
becomes thinner, leading to a larger thermal diffusivity just be- 
low it. This in turn implies stronger damping (see Eq.|6jl and the 



disappearance of low frequency waves. Furthermore, as stellar 
mass rises up to about 1.35 M & (i.e., T e ff ~6400 K), so does 
the flux associated with a given frequency. This is related to 
the rise of the luminosity with mass, and thus, of the energy in 
convective motions. For even more massive stars, the flux asso- 
ciated with a given frequency remains about constant, but most 
low frequency waves are damped too rapidly to exist. 

Let us recall that these results are obtained for stellar struc- 
tures on the ZAMS, and that we have found no significant de- 
pendence of the wave spectrum with age on the main sequence; 
it can thus be considered as constant for the duration of the 
main sequence. 

In this description, waves with frequencies up to N c (the 
Brunt- Vaisala frequency at the base of the convective envelope) 
can be excited. However, the high-frequency waves will have 
little impact on momentum evolution; indeed, as frequencies 
increase, filtering by the shear layer (see §3) becomes less effi- 
cient and less differential. Furthermore, damping is so small 
that it leads to only a small amount of momentum redistri- 
bution. Finally, as wave excitation diminishes with frequency, 
they carry much less momentum overall. These high-frequency 
waves will form standing waves, the g-modes of helio- and 
astero-seismology. 

2.3. Excitation by convective penetration 

Greater uncertainties exist on wave generation via convective 
penetration. The only theoretical estimates that exist have been 
made by Press (1981) and Garcia Lopez & Spruit (1991), and 
later used by Zahn et al. (1997). It is the formulation we shall 
adopt here. 

This description is based on the matching of wave pressure 
fluctuations with those of turbulent convection. Furthermore, 
it takes into account the combination of turbulent eddies of a 
given size to produce larger fluctuations. The range of horizon- 
tal scales available is thus 

/ co \ 3/2 

Q<£<( u with 4 = e cl—J (3) 
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Fig. 2. Momentum spectrum in £ and v for various masses for (left) Goldreich, Murray, & Kumar model and (right) Garcia Lopez 
& Spruit model 



where l c is the spherical harmonic number associated with the 
largest convective scale and u> c is the corresponding convec- 
tive frequency. Turbulence is assumed to follow a Kolmogorov 
spectrum, all frequencies with a> > a> c thus being available. 
The associated energy flux Te F is 

^ /( ^)^ _!__(!__) (4) 

and the corresponding angular momentum flux is still given by 
Eq. ©. 

2.4. A word of caution 

Let us stress that the fluxes calculated here are somewhat un- 
certain, especially in the second case. Indeed, it is well known 
that the "structure" of convection, which contains e.g. plumes 
that travel down across the whole convection zone, is not well 
reproduced by the mixing length model even though the con- 
vective velocities are more or less correct (e.g. Hurlburt et al. 
1986). However, their differential properties should have the 
proper dependence, and this will be our focus in the rest of this 
paper. 



3. Formation of a shear layer 

3.1. Differential filtering of the waves by a shear layer 

As described by several authors (Gough & Mclntyre 1998, 
Ringot 1998, KTZ99, Kim & MacGregor 2001), the most ob- 
vious feature of momentum transport by gravity waves is the 
formation of a double peaked shear layer just below the con- 
vection zone (see Fig. [3}. 

Gough & Mclntyre (1998) and Ringot (1998) argued that 
such a layer would prevent waves from propagating beyond. 
Indeed, local radiative and viscous damping is greatly depen- 
dent on the local frequency and increases when that frequency 
diminishes (this property leads to the formation of the dou- 
ble peaked shear layer). The local momentum luminosity in- 
tegrated over the whole spectrum writes 

■CAr) = J] ex P £)] . (5) 

The local amplitude exp [— r(r, <x, f)] depends on the integrated 
damping due to thermal diffusion Kj and (turbulent) viscosity 
V t 

T(r,cr,e)= (6) 
3 r r ° NNl I N 2 V dr 
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Fig. 3. Oscillating shear layer below the surface convection 
zone of a 1.2 M star. The dotted line shows the initial ro- 
tation profile. With the surface rotating slower than the core, 
a prograde shear layer is initially formed, followed by a ret- 
rograde one. When the shear becomes too intense, turbulent 
viscosity acts to merge the prograde layer with the convection 
zone, leaving behind the retrograde layer. A new prograde layer 
forms behind, and migrates towards the convection zone when 
the retrograde layer is absorbed. The cycle then resumes. 

where N 2 = Nj + N 2 is the Brunt-Vaisala' frequency, Nj is 
its thermal part and TV 2 is due to the mean molecular weight 
stratification (Zahn et al. 1997). 

The prograde peak 2 thus filters prograde waves, while the 
retrograde peak filters retrograde waves. If the peaks were in- 
finite in height, not a single wave could travel through it. This 
is however not the case, as the magnitude of the shear layer is 
self regulated by shear turbulence 3 . 

As described by TKZ02, if no initial differential rotation 
is present, the average wave momentum luminosity that tra- 
verses the shear layer is null, and there is no net effect on the 
interior. However, if differential rotation is initially present, as 
in the case of low mass stars that are braked by a magnetic 
torque, the average magnitude of the two peaks is not equal; in 
the case of interest (i.e. with the surface convection zone rotat- 
ing more slowly than the interior), the prograde peak is always 
larger, leading to differential filtering that favors the penetra- 
tion of retrograde waves. These waves are then damped in the 
whole radiative region that they spin down by depositing nega- 
tive momentum. 

3.2. Dynamics of the shear layer 

To study the effect of gravity waves in low mass stars, one 
must thus be able not only to evaluate the global momentum 
luminosity in waves, but also the dynamics of the shear layer 



2 The prograde peak corresponds to the layer that rotates more 
rapidly than the convection zone, while the retrograde peak designates 
the layer that rotates more slowly. 

3 Baroclinic instabilities could also set in, and would have the same 
effect. 




Fig. 4. Net luminosity vs time (in years) at 0.03 R* below the 
surface convection zone. Solid lines correspond to the case of 
strong initial differential rotation (6Q. = 0.05 yuHz over 0.05 R*) 
and dashed lines to the case of small initial differential rota- 
tion (<50 = 0.00005 //Hz over 0.05 R t ). Bold lines represent the 
mean luminosity, thin lines represent the instantaneous lumi- 
nosity. The curves correspond to a 1.2 M Q model. The shear 
layer oscillation occurs here with a period of ~ 1 year. 

and the resulting differential filtering. This is addressed here by 
performing numerical simulations of angular momentum evo- 
lution. We study a few cycles and evaluate the resulting net mo- 
mentum luminosity for various masses and various differential 
rotations. 

Angular momentum evolves under the action of gravity 
waves and shear turbulence according to 



d r 2 n Id 



4 <9Q 
pv,r — 
Or 



3 1 8 r ^ 
07T r L Or 



(7) 



where p is the density, Q the angular velocity and v, the tur- 
bulent viscosity. The local momentum luminosity -Cj(r) is in- 
tegrated over the whole wave spectrum. Equation (7) is solved 
throughout the radiation zone. The upper boundary condition 
expresses conservation of momentum of the star as a whole, 
with the convection zone rotating as a solid body. On their way 
to the core, waves are reflected when their local frequency cr 
equals the Brunt-Vaisala frequency 



cr(r, m) = u) - m [Q(r) - Q c ] = N. 



(8) 



Waves then deposit more momentum as they travel back to the 
surface convection zone. The high-frequency waves of lowest 
degree considered in this study can be reflected up to 5000 
times before they are completely damped; these contribute little 
to momentum redistribution 4 . 

Figure 3 shows the oscillation cycle of the shear layer for 
our 1 .2 M star. Each model (i.e., each stellar mass) has a differ- 
ent characteristic time-scale for this oscillation which depends 
mainly on the net wave flux in the high degree, low frequency 
wave. It is thus shortest for masses around 1.2 



1.3 M Q . 5 



4 As the frequency rises, damping diminishes and standing waves 
may form. 

5 Let us note that the very different time-scale found in TKZ02 for 
the Sun is related not only to the smaller flux but also to the reduced 
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Fig. 5. Net average (bold) and instantaneous (thin) luminosities 
0.03 R* below the surface convection zone after a few cycles 
for various masses and in the case of "strong" initial differen- 
tial rotation (6Q. = 0.05 fiHz over 0.05 R*). The dashed line 
corresponds to positives luminosities. 



4. Momentum extraction 

4.1. Fluxes 

In all our stellar models, we follow the detailed time evolu- 
tion of the shear layer according to Eq. (0 over several cycles, 
for two values of the initial differential rotation (a very small 
differential rotation of 8£l = 0.00005 juHz over 0.05 R», and a 
larger differential rotation 5Q. = 0.05 fjHz over 0.05 R*). We 
calculate the net luminosity below the filtering shear layer, at 
r = R cz - 0.03/?*. Resulting instantaneous and mean luminosi- 
ties are detailed for the 1.2 M model in Fig. [4] The change 
of sign of the instantaneous luminosity illustrates the varying 
asymmetry of the double peak shear layer. In the case of very 
small initial differential rotation, the curve is symmetric and the 
mean luminosity is always close to zero. On the other hand, if 
the initial differential rotation is large, retrograde waves are fa- 
vored, as illustrated by the negative mean luminosity. It is this 
net mean luminosity below the shear layer that gives rise to a 
net momentum extraction in the interior. 

In order to obtain differential filtering, it is important to in- 
clude several different modes; 



number of frequencies and degrees used in order to be able to perform 
the solar calculation over the required 10 7 time-steps. 



^Efficient transport 
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Fig. 6. Net luminosity at 0.03 R* below the surface convection 
zone as a function of r e (j(zams) for an initial differential ro- 
tation of 6Q. = 0.05 yuHz over 0.05 R*. Triangles: without the 
Coriolis force; circles: including the Coriolis force. The verti- 
cal lines correspond to the requirement of angular momentum 
efficiency as described by TC98 (they correspond to the hori- 
zontal lines of their Fig. 2). 



- High-£ modes, that damp close to the convection zone (see 

- Low-£ modes, that penetrate deeply. 

Furthermore, it is the lowest frequency waves (which undergo 
the largest amount of differential damping) that give rise to mo- 
mentum extraction. While wave frequencies range from ~ to 
N, for each considered mass, there is a maximum frequency for 
waves that contribute to differential filtering. That maximum 
frequency has been determined for each individual mass (see 
Table 1). In the most massive models considered here, which 
correspond to stars inside the dip, this maximum frequency is 
much higher than in lower mass stars. Indeed, as can be seen 
from the wave spectrum (Fig.|2ji, in those more massive mod- 
els, high-^ modes of low frequency experience too much damp- 
ing to travel beyond a single vertical wavelength in the vertical. 
One must thus include higher wave frequencies in order to cap- 
ture high-^ modes that produce the oscillating double peaked 
shear layer. 

The Garcia Lopez & Spruit model for wave excitation by 
convective penetration leads to the generation of larger fre- 
quency waves, leading to a limited amount of extraction. This 
is especially true for the more massive stars presented in this 
study and is thus less relevant to our discussion here 6 . For that 
reason, that model has been ignored in detailed calculations. 

The time evolution of the resulting filtering in all stellar 
masses considered is shown in Fig. [5] The abscissa represents 
the time iteration, each mass having its own characteristic time- 
scale for oscillation. The results are thus shown in terms of this 
characteristic time, and not in a real physical time sequence. 
The net luminosity slightly increases with the stellar mass up to 
1 .2 M , and then dramatically decreases as one moves to stellar 
masses inside the Li dip. For the largest mass considered here 



Let us remind that both generating sources are actually additive. 
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Table 2. Effect of the Coriolis force. 
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(1.44 and 1.45 M ), the process of differential filtering does 
not lead to momentum extraction anymore (the net luminosity 
is actually positive although very small for these stars on the 
blue side of the dip 7 ). 

Figure [6] represents the average luminosity after 500 time 
steps as a function of effective temperature. On the same figure 
are shown the limits established by TC98 in terms of momen- 
tum extraction required in order to consistently understand the 
Li dip in terms of rotational mixing. 

4.2. Coriolis force 

In order to fully describe wave transport in rotating stars, we 
must also take into account the influence of the Coriolis force 
on waves. Its first order effect is to change the horizontal struc- 
ture of the modes, confining them closer to the equator. The 
maximum co-latitude of propagation # cr jt is given by the condi- 
tion 

a- 2 = 4Q 2 cos 2 6U, (9) 

(KTZ99). For a given rotating velocity, this condition implies 
that only a fraction 

S= sin6»d0 (10) 

of the surface will support waves and the corresponding lu- 
minosity should be diminished accordingly. Table [2 gives the 
critical angle for propagation and the corresponding efficiency 
for a typical frequency of 1 fiHz. The rotation velocities used 
here correspond, for each stellar mass, to the mean velocity ob- 
served at the age of the Hyades. Let us note that the fraction of 
the surface S varies linearly with v, and that the whole surface 
will support waves when v > 2 v rot . 

Exact calculations can be made by reducing the wave lu- 
minosity corresponding to each frequency by the proper coef- 
ficient. This has been done for 3 models (1.2, 1.3 and 1.4 M Q ). 
The amplitude of low frequency waves (that dominate the for- 
mation of the shear layer) being reduced, including the Coriolis 

7 This reverse effect is probably linked to the fact that thermal 
damping dominates just below the convection zone. Then, the low-^, 
high frequency waves that are included and that travel back and forth 
over 100 times, are damped mostly there. However, retrograde waves 
that have a larger frequency in the interior are less damped there than 
the prograde waves, leading to an effect opposite to what is observed 
in other stars. 



force leads to a somewhat smaller double peaked shear layer, 
that oscillates on a slightly longer period. Furthermore, dif- 
ferential filtering corresponding to momentum extraction is 
also dominated by rather low frequency waves. This leads to 
a smaller efficiency of momentum extraction by waves (see 
Fig.|6j. However, the overall mass dependence is not affected. 

5. Discussion 

Several abundance anomalies are best explained in terms of 
large-scale mixing. As a natural physical parameter (beyond 
mass, chemical composition and age), rotation is an important 
factor in understanding various features which are not predicted 
from classical stellar models (see Maeder & Meynet 2000 for 
details). Furthermore, if one is to built a consistent model of 
rotational mixing, it has to apply to all stars. 

Recent improvements of the description of the Eddington- 
Sweet meridional circulation by Zahn (1992) and Maeder & 
Zahn (1998) have been applied to several kinds of stars. These 
studies show that, with the same free parameters, it is possi- 
ble to explain abundance anomalies in B (Talon et al. 1997) 
and O stars (Maeder & Meynet 2000) as well as the blue side 
of the Li dip (TC98, PTCF03) and Li anomalies in sub-giants 
(Charbonnel & Talon 1999). 

However, in order to understand the helioseismic rotation 
profile, some extra mechanism for momentum transport is re- 
quired in low mass stars (Matias & Zahn 1998). TC98 argued 
that, for rotational mixing to be responsible for the Li dip, this 
mechanism has to become efficient at an effective temperature 
around ~ 6500 - 6600 K, i.e., in a region where the stellar con- 
vective envelope becomes substantial. If the same description 
is to be applied to all stars, there must be some physical reason 
for the lack of efficiency of this mechanism in hotter stars. The 
goal of this paper was to establish whether gravity waves do 
have the proper effective temperature dependence. 

The calculations we performed here show that the built up 
of a double peaked shear layer allowing differential filtering of 
the waves is possible in low-mass main sequence stars up to an 
effective temperature of ~ 6700 K. The net momentum lumi- 
nosity below the shear layer increases with stellar mass up to 
1 .2 M (i.e., r eff ~ 5900 - 5950 K); it then drops in hotter stars, 
even becoming slightly positive for stars more massive than 
~ 1.44 M (i.e., r eff > 6700 K). While the quantitative results 
shown here are based on the excitation model by Goldreich et 
al. (1994), the Garcia Lopez & Spruit (1991) model would lead 
to an even stronger decrease in efficiency. Indeed, as can be 
seen in Fig. 12 in the more massive models, there are simply no 
low frequency waves, implying an even more abrupt decrease 
in the efficiency of differential filtering. 

Let us discuss a final issue. As can be seen from Fig. 1, 
the dispersion in lithium abundances is somewhat larger on the 
cold side of the dip (T^g ~ 6000 - 6500 K) than in cooler stars. 
A very straightworward explanation of this arises within our 
framework. As we have shown indeed, for a given stellar mass 
a higher rotation velocity implies a stronger reduction of the net 
momentum luminosity in waves due to the Coriolis force. As 
a consequence, the efficiency of angular momentum extraction 
by the waves is somewhat lower within faster rotators, leading 
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to a larger efficiency of the hydrodynamical processes, and thus 
to stronger Li destruction. In this context, the larger lithium dis- 
persion measured in stars laying immediatly on the right of the 
dip can be interpreted as a consequence of the larger dispersion 
in rotational velocities that these stars do exhibit in comparison 
with the less massive ones (Gaige 1993). Detailed calculations 
would have to be performed to check whether this effect could 
be large enough to explain the observed dispersion. 

All our results clearly indicate that momentum trans- 
port by gravity waves has the proper mass dependence to 
be the required process in low-mass main sequence stars on 
the cold side of the Li dip. 

Note that a similar mass dependence is not expected if mo- 
mentum transport is dominated by magnetic fields. There is 
another test which could permit to differentiate the magnitude 
of transport by magnetic field compared to transport by waves 
which applies to the Sun. Certain early helioseismic inversions 
(Elsworth et al. 1995, Corbard et al. 1997) suggested that the 
central region of the Sun (below ~ 0.4 R Q ) could rotate some- 
what slower than the rest of the radiative region. These results 
seem to be confirmed by recent SOHO data (Couvidat et al. 
2003). TKZ02 explained how such a feature could be produced 
by gravity waves. Both arguments favor momentum transport 
by gravity waves with respect to magnetic field. 

We are now in a position to present a coherent picture of 
rotational mixing in main sequence stars of all masses: 
On the hot side of the Li dip and in more massive stars, the 
transport of angular momentum and of chemicals by merid- 
ional circulation and shear instabilities do explain the Li as 
well as the He and CNO patterns. In lower mass stars, gravity 
waves dominate the transport of angular momentum, thereby 
reducing the magnitude of meridional circulation and shears 
and shaping the Li pattern on the cold side of the dip. Within 
this framework, we thus predict that Pop I main sequence stars 
with initial masses lower than ~ 1.4 M Q must be quasi-solid 
body rotators, as the Sun is. 
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